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Signatures of critical behaviour in nuclear fragmentation are often 
based on arguments from percolation theory. We demonstrate with 
general thermodynamic considerations and studies of the Ising model 
^ ' that the reliance on percolation as a reference model bears the risk of 

missing parts of the essential physics. 



^ ' 1 Introduction 



^ . The existence of phase transitions in nuclear matter was proposed more than 20 

years ago in conjunction with the description of the structure of neutron starsS'i. 
This led to the question whether phase separation!, dynamical instabilities or 
critical phenomenal play a role in the disintegration of highly excited nuclei pro- 
duced in proton or heavy ion induced collisions. The observation of a power law 
for the fragment size distribution i'i was the first experimental hint for a second 
order phase transition, but the determination of critical exponents! or other pa- 
rameters! turned out to be difficult with the mostly inclusive data available at 
that time 000. In the past 3 years, however, data from a new generation of 
experiments with an almost complete coverage for the decay products of an ex- 
cited system allowed for rapid progress. The observation of potential signals of a 
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first order phase transition EJo and the attempt of a quantitative determination 
of critical exponents sparked an intense discussion on what rehable signatures 
for a first or second order phase transition in a small and dynamically evolving 
system are. 

In the following we will take the recent claim lillii, that a continuous phase 
transition with critical exponents consistent with those of a liquid-gas system 
has been observed in the fragmentation of nuclei, at face value and revisit the 
theoretical arguments used in this analysis. 

2 Critical Phenomena in a Nutshell 

Near a critical point, the order parameter cj), associated susceptibility specific 
heat c, and the correlation length ^ exhibit a power law dependence on two 
reduced control parameters e and h: 

0{e,h=o) = Be^ Liquid-Gas Magnet Percolation 

0(,=o,h) ~ h^'^ (j> p-pc M Poo 

Xi.,H=o) = re-- with X H Es'n, 

C{e,h=o) = ^e~" e ^ ^ pc-p 

^{e,h=o) = ^oe"'' h p- po{T) H Pg 

For liquid-gas, magnetic and percolation systems the order parameter is the 
difference between liquid and critical density, magnetization per spin M or the 
fraction of sites in the largest cluster Pqo, respectively. The susceptibility x corre- 
sponds to the isothermal compressibility kt = and the second moment of 
the cluster size distribution for the liquid-gas and percolation case. The 'thermal' 
control parameter e is the difference of the temperature T or bond probability p 
to the critical value while the 'field' control parameter h is given by the chemical 
potential p, the external field H or the ghost bond 00 probability Pg. 

3 Prom Thermodynamics to Percolation and Back 

The exponents (3 and 7 can be determined from a measurement of (/>(e,/i=o) and 
X{t,h=o)- For liquids it is indeed possible to determine Ap = p — Pc and |^ directly, 
e.g. by measuring the density profile in the gravitational field 000. 

In the nuclear physics case, the only directly measured quantities are the size 
and momenta of clusters produced in an interaction. The standard argument 
to deduce Ap and kt from the cluster size distribution uses the Fisher droplet 
mo deli. In this Ansatz the grandcanonical partition function Z is expressed as 
a sum over cluster yields Y(a) which depend on a surface energy s and chemical 
potential difference between liquid and gas phase p: 

lnZ = Y^ Y^A) with = qoA-^ exp (-i {sA'' + pA)) 
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From this follows immediately that kt is proportional to the second moment of 
the cluster size distribution 

It can also be shown il that Ap is proportional to the fraction of constituents 
belonging to the largest cluster Poo- This strong similarity to percolation, where 
Poo and the second moment also play the role of order parameter and susceptibil- 
ity, has led to a widespread use of percolation in the modeling and interpretation 
of nuclear physics experiments. In particular, percolation has served as reference 
model for the development of methods to extract critical exponents 

Even though there is a mathematical connection between percolation and 
the thermodynamics of interacting systems (the g — 1 limit of the Potts model 
corresponds to bond percolationE3) there are also some significant differences: 

• Percolation is usually discussed in terms of only one control parameter, the 
bond or site probability, leaving the impression that the 'field' control param- 
eter, corresponding to chemical potential or average density in interacting 
systems, is of minor importance. 

• Finite size effects depend strongly on the ensemble in interacting systemsilSil, 
but there is no direct equivalent to the concept of an ensemble in percolation. 

• Percolation theory simply starts with the definition of what a cluster is. Ther- 
mal systems, on the other hand, are usually defined in terms of the interaction 
between their constituents and the appropriate definition of a cluster is, as 
will be shown, nontrivial. 

The almost exclusive reliance on percolation as a reference model bears therefore 
the risk that important parts of the physics are missed. In the following we will 
consequently turn back to the thermodynamic basics and use the Ising model as 
guidance. 

4 Prom Constituents to Clusters 

As a first step it is interesting to establish the connection between kt and the 
cluster size distribution with minimal assumptions and without using a specific 
model. We follow an idea given by Alexandrowicz il for the Ising model but 
generalize it to an arbitrary system. In the grand canonical ensemble, kt is 
related to the fluctuation of the particle number cr^ by 



2 l.T^^ 



= kT — KT 
T,V ^ 



In the following, we consider a system with Ns constituents coupled to a reservoir 
with Nji constituents: 
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. gygi^gj^ / ^ Define for constituent i: 
A^: Reservoir | I Ns I | _ J 1 if in 5 
Nr ' ' J ~ 1 if in i? 



To express in terms of clusters we assume that the constituents are grouped 
into clusters with the properties: 

PI A cluster is either completely in S or in R. 

P2 Constituents in different clusters are uncorrelated. 

can now be rewritten as a g correlation. The sum over constituent pairs can 
be split into two parts, one where the pair is in the same cluster and one where 
it is in different clusters: 



^2 





{Nsf 






\ ij 




\ ij 





with q — -^^^ 



N 



diff 



(P2) 

= ' 



S R 

In the last line, the sum over constituents was rewritten as a sum over clusters 
E l*^*!^ ^P^^^ ^'^^'^ parts, running over S and i?, respectively. In 

ij Ci 

same 

the thermodynamic limit, Nji ^ oo , g ^ 0, the second term goes to zero if the 
largest cluster in R grows slower than Nr (in other words: No condensation in 
R) and one finally gets 

kT {Nsf kTp^ ' ^ ^ 

where Ug is the concentration of clusters of size s. This is equivalent to Eqn. (|l|) 
but was derived with only one essential assumption, that constituents in different 
clusters are uncorrelated. 



5 The Ising Model 

One of the simplest thermodynamic models with a phase transition is the Ising 
model, given by the Hamiltonian: 

n = -J (r,a,-Hj2cr, (3) 

<i,j> i 
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It can be interpreted as a model for a magnet, with a = ±1 representing up 
and down spins, or as a lattice model of a gasEO where a = ±1 now indicates 
whether a site is occupied or empty. The number of sites A^, magnetization M 
and external field H correspond in the lattice gas interpretation to the volume 
V , density p and chemical potential /x, respectively. 

Order parameter (|M|) and susceptibilities x aiid x' fo^' T > and T < T^, 
respectively, are defined in terms of constituents as: 









X = 


f (M^) = 




x' = 




-{\M\f) 



(4) 



M and x can be expressed in terms of clusters using arguments which are for x 
analogous to the ones in the previous section. A detailed derivation along those 
lines was given by de Meolll while other authors arrived at the same result from 
different viewpoints 00110 One obtains for the cluster observables M, x a^id 

x'- 

^ (^oo) < (|M|) 



M 
X 
X' 



kT 

^ (x; + N{Apy 



> 



X 

x' 



(5) 



with the fluctuation of largest cluster (AP)^ 



m 



< rim > where 



(P^) — (-Poo)^ and the percolation 
is the number of clusters per site 



susceptibility Xp = 
with magnetization m. 

While the order parameter is indeed given by the relative size of the largest 
cluster, we observe two significant differences between x and Xp- 

• The trivial ^ factor (see also Eqn. and (0)) makes sure that the paramag- 
netic or ideal gas limit x ^ is approached for T ^ oo where Xp 1- Even 
though it does not affect the asymptotic power law behaviour for T Tc it 
will change the effective exponents if a wide temperature range is considered, 
like in the nuclear exponent analysis 0. 

• The susceptibitity x' for T < T^is given by the sum of x'p, the second moment 
taken without the largest cluster, and the fluctuation of the largest cluster 
(AP)2. 

Numerical results for the 3d Ising model with periodic boundary conditions are 
shown in Fig. |l|. The calculations were done with the Swendsen-Wang and 
Wolffil algorithms in implementations similar to Ref.EsI using a random number 
generator proposed by Ziff 00. In the left two columns we compare the cluster 
observables to the true values given by Eqn. (^. Consistent with Eqn. (|^) we 
find M slightly smaller than (|M|) and x equal to x while x' is only close to x' 
for T < Tr. 
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Figure 1: Left two columns: Comparison of cluster observables (Eqn. (|^), full dots) and true 
values (Eqn. (^), full lines) for anlaing model with 8'^ sites. The dotted lines indicate the bulk 
critical temperature Tc — 4.51152E3. 

Right column: Susceptibility ratio x/x' *t Tc as a function of the linear system size L with 
(full dots) and without (open dots) the (AP)^ contribution to x'- 



It is interesting to note, that the main contribution to x' is the fluctuation 
term (AP)^ while the second moment term -pf^x'p (indicated by open dots in the 
top middle frame of Fig. ||) carries only a small part of the signal. This can be 
seen more quantitatively in the right panel of Fig. |1] where the susceptibility ratio 
x/x' at Tc is shown as a function of the linear system size L with (full dots) and 
without (open dots) the (AP)^ contribution to x'- Depending on system size, 
between 60 and 75% of the susceptibility are carried by the fluctuation term. 



6 Critical Clusters 

Clusters can be defined in interacting systems in many different ways. We call 
the result of a definition 'critical clusters' if order parameter and susceptibility 
correspond to size of largest cluster and second moment, or formally, if Eqn. (|^) 
holds. We saw in the previous two sections that a necessary condition is the 
noncorrelation of constituents in different clusters but have not given so far an 
explicit definition. 

The most obvious way to define a cluster is the 'geometrical cluster' which 
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corresponds to a domain in the Ising model, or in the lattice gas picture to a 
connected region of occupied sites surrounded by empty sites. However, it is 
easy to see that this violates the noncorrelation requirement and it has indeed 
been shown that in 3d 'geometrical clusters' lead to a percolation transition at 
temperatures below the thermal transitioned, thus do not fulfill Eqn. (|^). 

The proper prescription for 'critical clusters' for the Ising model was given by 
Coniglio and Klein e3 and requires that the geometrical clusters are broken into 
smaller pieces with a bond percolation with pb = 1 — exp(— |^). 

The definition of a 'physical cluster', based on pairwise binding, proposed 
by HillS and recently used in a lattice gas model for nuclear fragmentation @, 
results in bond probabilities quite similar to the Coniglio-Klein values. Even 
though both definitions give numerically similar results for small system sizes 
they are clearly not equivalent in the thermodynamic limit. 

7 Summary 

Percolation describes remarkably well the cluster distributions and correlations 
in nuclear fragment at ion 00. This, together with its simplicity and ease of use, 
made it the ideal reference model for the study of signatures of critical behaviour. 
But given the hints, that we observe a liquid-gas rather than a percolation phase 
transition, we have to face the limitations of percolation. Even scratching at the 
surface of a proper thermodynamic description raises many issues: 

• Critical clusters have rather remarkable properties, they are neither well 
separated nor compact ||, are interacting but have nevertheless no correlations 
between constituents of different clusters. So one might wonder whether crit- 
ical clusters are mere mathematical constructs, like the clusters in Mayer's 
cluster expansion, or real physical entities. All attempts to infer signatures of 
critical behaviour in nuclear physics not only imply that critical clusters are 
observable objects, they also assume that the distribution of clusters formed 
in the decay of a system is representative of the equilibrium distribution at 
some freeze-out condition. Even though there are attempts to support this 
connection with model calculations e3 it remains to be seen whether there is 
a more rigorous way to justify this. 

• Ensembles: Although the experimental situation is certainly better repre- 
sented by a microcanonical treatment we used for simplicity for all arguments 
in sections ^ and ^ the grand canonical ensemble. While the expectation values 
of extensive quantities do not depend on the ensemble in the thermodynamic 
limit it is easy to show that this is not the case for fluctuationsi^. The cluster 
observables (Eqn. (|^)) are therefore likely to be ensemble dependent. 

• Role of control parameters: The exponent analysis rests on the precon- 
dition that the temperature of the system is varied while the second control 

^they are in general fractals with a dimension yh = ^ or about 2.5 in the 3d Ising class 



7 



parameter, the density, stays constant and close to its critical value. Assuming 
that we observe conditions along a freeze-out line and taking Papp's schematic 
model0'i guide one sees that the reverse might be true, that the temper- 
ature is almost constant but that the density varies. In this case one would 
determine different exponents, e.g. 1/6 rather than (5 (see section |^). One 
should also keep in mind that the relative size of the largest cluster correlates 
with the density of the liquid phase (thus the order parameter) only on the 
critical isochore but is in general more readily interpreted as a measure of the 
mass fraction of the liquid phase. 

• Corrections to scaling: The power laws listed in section |^ hold only asymp- 
totically for e ^ while a description in a wider range is possible with correc- 
tion to scaling termsii or a crossover approach 3. The size of the critical region 
depends strongly on the exact form of the interaction and is rather small in 
real liquids but substantially larger in models with only next-neighbor or 
short range interactions. The leading correction term for the order parameter 
in the Ising model S is not only substantially smaller in magnitude but has 
even a different sign as compared to typical liquids^. It is therefore uncertain 
whether simple schematic models, like percolation or Ising, are adequate for 
the modeling of finite size and finite control parameter effects in the nuclear 
case. 

• Field gradients: The hallmark of critical phenomena is scale invariance 
which at least requires homogeneous conditions throughout the system un- 
der study. A Coulomb or a radial flow ^ field will prevent the growth of 
fluctuations across the whole system much like the chemical potential gra- 
dient caused by the gravitational field is limiting the usable sample size in 
earthbound experiments on liquidsi§. 

All those points will have to be addressed before the experimental results can be 
connected to the parameters of bulk nuclear matter in a quantitative way. 
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